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Abstract 

We have numerically investigated the mean-field dynamics of the the p-spin 
interaction spin glass model with p=3 using an efficient method of integrating 
the dynamic equations. We find a new time scale associated with the onset of 
the breakdown of the fluctuation-dissipation theorem in the intermediate time 
regime. We also find that the off-equilibrium relaxation exhibits a sub-aging 
behavior in the intermediate times and crosses over to a simple aging in the 
asymptotic regime. 
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Since the recognition [||] of possible deep connection between structural glass and some 
class of spin glass systems such as the p-spin model fHH, the spherical p-spin model [|J, 
due to its analytic accessibility, has been the subject of intense research in both statics and 
dynamics. Here we examine some open questions (such as scaling behavior) || of the off- 
equilibrium dynamics of the mean field spherical p-spin model by developing a new method 
of integrating the dynamic equations of the model. 

We consider the system of N spins Si, S2, ■ ■ ■ Sn which interact via Hamiltonian 

H — — 2_j Ji 1 i 2 ---i p Si 1 Si 2 ■ ■ ■ Si p (1) 

l<ii<j 2 --<i P <iV 

where the spins are continuous variables subject to the spherical constraint J2iLi Sf = N. 
The coupling Ji^...^ is a Gaussian random variable with zero mean and variance [Ji li2 ... ip ] j = 
p\/2N p ~ 1 . It is well known M that this model (with p > 2) exhibits an equilibrium phase 
transition at a finite temperature T = Tg from paramagnetic phase to spin glass phase 
characterized by one-step replica-symmetry breaking. This temperature is lower than the 
than the dynamical freezing temperature To, below which the dynamics is non ergodic. 

For the time evolution of the system, we consider the following dissipative Langevin 
dynamics 

T-%S t (t) = -z(t)S t (t) - Q^j + rk® (2) 

where z(t) is a Lagrange multiplier to enforce the spherical constraint for all times. In order 
to satisfy the detailed balance, the thermal noise f]i(t) is chosen to be Gaussian with zero 
mean and variance {Vi{t)Vj{t')) = 2r _1 T5jj(5(t — t'). The inverse of the kinetic coefficient T 
sets the microscopic time scale, which henceforth is set to unity. 

This Langevin dynamics governs the time evolution of the system starting from, for ex- 
ample, the random disordered spin configuration toward thermal equilibrium. Two physical 
quantities of interest, which quantitatively characterize this time evolution, are the two-time 
correlation function C(t,t w ) = Y,i=i[< Si(t)Si(t w ) >}j/N where (• • •) and [• • -]j represent 
averages over thermal noise and random coupling, respectively, and the response function 
R(t,t w ) = J2fLi[ ^h.fp)" where hi{t w ) is the external field turned on at time t w . In 
equilibrium, these quantities become time-translation invariant, i.e., C(t,t w ) = C(r) and 
R(t,t w ) = R(t), t = t — t w . Moreover, there exists a fundamental relationship between the 
correlation and the response, known as the fluctuation-dissipation theorem (FDT), which 
takes the form for the present irreversible dynamics 

R(T) = JM?m (3) 

where 9{r) is the unit step function which reflects the causality for the response function. 
We will see below that the non-equilibrium dynamics of the present system manifests an 
explicit strong dependence of the two times t and t w (aging) and an interesting modification 
of FDT . 

At this stage, we find it very useful to introduce an integrated response function F(t, t w ) 
defined as 
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F(t,t w ) = - f dsR(t,s). (4) 

J t w 

From (D F(t,t~) = and dF{t,t w )/dt w = R(t,t w ). FDT (g) then takes the form F(r) = 
9(r)[C(r)-l]/T. 

In the limit of N going to infinity, one can treat the dynamics exactly using the standard 
functional method such as Martin-Siggia-Rose formalism [|X0| Jlf| . In particular, the dynamics 
is governed by the closed set of coupled equations of C and F: 

d t C(t, t w ) = -z(t)C(t, t w ) + P{p ~ l > f dsC p - 2 (t, s)(d.F(t, s))C(t w , s) 

I Jo 

+ l f w dsC p - 1 (t,s)d s F(t w ,s), (5) 
2 Jo 

d t F{t,t w ) = -1 - z{t)F{t,t w ) + P{P ~ l) f dsC p - 2 (t,s)(d s F(t,s))F(s,t w ). (6) 

z Jt w 

where z(t) = T — pE(t), E(t) being the average energy per spin. The energy density 
E{t) is related to C and F as E(t) = -(p/2) J * ds C p ~ l {t, s)d s F(t, s). The same type of the 
equations have been derived in various physical contexts such as the dynamics of a long-range 
superconducting wire network (p = 4) [TT|, the dynamics of Amit-Roginsky model(p = 3) 
fj~2|, and the dynamics of a particle in the random potential in large dimensional space 



T3HI6|. 



As Cugliandolo and Kurchan (CK) || have shown, the dynamics exhibits the two distinct 
regimes depending on the relative magnitude of the two times r(t = t + t w ) and t w . The 
first regime (while t,t w — * oo, the time difference r = t — t w is finite, i.e., r/t w — > 0) is the 
regime where the time-translation invariance and FDT hold. That is, C(r + t w , t w ) = C(r), 
F{t + t w ,t w ) = F(t) and TF{r) = (C(r) — 1) independent of a given t w . In this regime, 
the dynamic equation for C(r) can be easily derived from the equation for the integrated 
response and is given by 

(d T + T)C(t) +p(e oo + (1 - C(t)) + ^ dsC p -\r - s)^ = 0. (7) 

where E^ is the long time limit of the energy density E(t), E^ = lim^oo E(t). One can 
recognize that this equation is quite similar to a schematic model developed in the context 
of the mode-coupling theory (MCT) for the glass transition in structural glasses [ 17| . In 
particular, apart from the term involving (1 — C(r)), (^) with p = 3 is identical to the 
Leutheusser model (without the inertial term) Jl8| . Hence one can sense that this equation, 
as in the Leutheusser model, may lead to a dynamic transition to a nonergodic phase where 
the long-time limit of C(r) is non- vanishing. It is easy to show from ([/]) that the non- 
ergodicity parameter q = lim,--^ C (r) is then related to -E^ via 

- (T - pE..) + ^(1 - if- 1 ) = (8) 

But the full understanding of the FDT dynamics requires the asymptotic value of the energy 
density, E^, for which one has to consider the dynamics in different regime (aging regime); 
the two regimes are closely coupled to each other. 
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When the two times t and t w are large and well separated, the relaxation is very slow, 
and hence the time derivatives in (|5|) and (§) can be ignored. In this situation, it is found 
that the scaling Ansatz for C and F 



C(t, t w ) — C 



h{t v 



leads to 



h{t) 
p{p — 1 



F(t,t v 



h(t v 



hit) 



2T 2 



(i- q y q 



2 nP~ 2 



(9) 



(10) 



TF{\) = xC(X) - [1 - (1 - x)q], (11) 

and 

^co = ~ [1 - (1 - a;)^ (12) 

with x = {p — 2){l — q)/q for non vanishing q. The equation flTT| ) with x < 1 is a modification 
of FDT in aging regime and is one of the most important results obtained from the asymp- 
totic analysis. Note that although the actual dynamics will select the form of h(t) uniquely, 
the above asymptotic analysis holds for an arbitrary monotonically increasing function h(t) 
(known as the time-reparametrization invariance): the function h(t) remains undetermined 
within the asymptotic analysis. 

We now go back to the equation (0) and discuss the dynamics of FDT regime. First 
note from (|8|) that q = is always a solution for all temperatures. We see from flTU| ) that 
the non- vanishing q starts to appear at the temperature T* = [2p 1 ~ p {p — l)(p — 2) p_2 ] 1 / 2 
at which q = q* = (p - 2)/p (q* = 1/3 and T* = 2/3 for p = 3). Thus below T\ starting 
from the initial state with zero energy density |19] {e.g., the state with all spins up), the 
system always chooses the solution with higher energy (which is the highest TAP state 
at a given T [p|,p0|-|22]] ) among these two solutions. Thus for T D < T < T*, q = is 
the genuine solution, and E^ = — 1/(2T), x = 1. The system is ergodic. The dynamic 
transition temperature is determined by x{T = Tp) = 1 for q ^ 0. This condition 
with (0) lead to q D = {p - 2)/{p - 1) and T D = [p{p - 2) p - 2 /{2{p - If- 1 )} 1 ' 2 (for p = 3, 



Td = y3/8 ~ 0.612 ■• • and q D = 1/2). Therefore the term involving (1 — C(r)) drops 
out in ([7p and the resulting equation is the same as a schematic mode-coupling equation 
for supercooled liquids. Therefore the relaxation dynamics above the dynamic transition 
exhibits the well-known scaling laws derived with MCT f23|. For T < Td the system 



chooses q ^ 0, x < 1, and E^ + 1/(2T) > 0. Hence the system is non-ergodic. Thus the 
term p{E 00 + 1/2T)(1 — C{r)) in (0) is turned on. This additional term makes the dynamics 
in the FDT regime in the present model differ from that of the MCT for supercooled liquids. 
In particular, near the transition q{T) shows a linear behavior q{T) = qo + const. {Td — T) 
instead of a square-root singularity observed in MCT. Also whereas the critical relaxation is 
seen very near and at the transition in MCT, in the present situation the correlation exhibits 
a critical relaxation C(r) = q + const. r~ a for all temperatures below the transition with the 
exponent a related to the FDT violating factor x as T 2 [l — a]/r[l — 2a] = x/2 15|| , T being 
the gamma function. 
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Though the asymptotic analysis put forward by CK provided quantitative informa- 
tions on the both FDT and the asymptotic aging dynamics, it leaves unexplored the non- 
equilibrium dynamics in the intermediate time regime. Moreover, the time reparametrization 
invariance leaves undetermined the form of the function h(t) which the original set of causal 
dynamic equations would select in the course of its evolution. Therefore, we do not know 
what type of scaling feature the long time aging dynamics might exhibit. 

In this work, we have developed an efficient way of integrating the mean-field dynamic 
equations (§) and (|6]), which can extend the solution to the time regime long enough so that 
one can clearly investigate the aforementioned important open questions. Basic idea is that 
since the relaxation becomes very slow at long times one can employ an adaptive integration 
time step fl24 |. There are two crucial ingredients in the present method. First, one has 
to work with the integrated response instead of the response function since the integrated 
response relax much more slowly than the response (which is basically a time derivative 
of the correlation). Another ingredient is to increase the integration time step along the 
both two time directions in such a way that the fast relaxing time regime should have many 
integration points. The details of the present method will be given elsewhere. 

Now we turn to the discussion of the numerical results. All presented results are for 
p = 3. Here we focus on the dynamics below the transition. Shown in Fig. 1(a) is C and F 
for T = 0.5. Both functions exhibit a strong t-dependence (aging effect) which persists up 
to longest time t. In contrast to the situation above the transition, here the aging will not 
be interrupted, and the system remains perennially out-of-equilibrium. The parametric plot 
—TF(t,t — t) versus C(t,t — r) for T = 0.5 shown in Fig. 2(b) exhibits several interesting 
features. After some time t ~ 10 2 , independent of time t, the FDT is established in the range 
q < C < 1. Then, after still longer time t ~ 10 5 , all the curves with different times t merge 
into a single straight line with slope close to the FDT violation parameter x ~ 0.565 • • •. 
Note the huge time interval (more than 9 decades!) it takes to reach the asymptotic regime. 
Therefore in the asymptotic regime, the FDT and its violation can be characterized by the 
two straight lines with slopes 1 and x meeting at C = q: —TF = 1 — C ior q < C < 1 and 
— TF = x(l — C) + (1 — x)(l — q) for C < q. This is in accordance with the result of the 
asymptotic analysis ([TT|). Similar breakdown of FDT has been observed in simulations of 
supercooled liquids |?5H 2~7| and spin glass [|8 . 



One interesting new feature we find in the off-equilibrium dynamics in the intermediate 
time region is the existence of a new time scale associated with the onset of FDT violation. 
Note that the breakdown of FDT occurs sharply at C = q (see Fig. 2(b)). Hence we are 
able to measure this time scale by reading off the times r p defined as C(t,t — r p ) = q for 
each fixed t. The double- log plot r p versus t in the inset of Fig. 2(a) indeed demonstrates 
that r p shows a power law dependence on t as r v ~ with ~ 0.68. Furthermore, the big 
difference in t p and t for large t implies that t p ~ with the same exponent. Since is 
smaller than 1, we find a new time scale characterizing the aging dynamics in intermediate 
time region. It is thus highly desirable to have a theoretical development which can provide 
a detailed information on the dynamics in the intermediate times and on the crossover to the 
asymptotic time regime. The presence of this time scale in the spherical SK model (p = 2) 



was demonstrated by Zippold et al [^9[ . 

We now discuss the scaling behavior in the aging regime. Recall that the asymptotic 
analysis loses information on the unique selection of hit). The scaling property in our 
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numerical solution is important since it will be able to determine the specific form of the 
function h(t) that the system actually selects. Figure 2(a) shows the relaxation of the 
correlation function C(r + t w ,t w ) for different waiting times t w at T — 0.5. Due to the 
adaptive integration steps along the both t and t w axes, it is inevitable that the short time 
data for each fixed t w is lost. It is natural to first examine the possibility of h(t) being 
a power of t, h(t) ~ t 1 ' . This means that the correlation in the aging regime shows a 
scaling behavior Ca{t + t w ,t w ) = C(r/t w ) (known as the simple aging). To check this 
scenario, we plot in Fig. 2(b) the correlation function C(r + t w ,t w ) against the rescaled 
time r/t w . Indeed, we observe that the relaxation data with longest waiting times collapse 
onto a single scaling curve. Thus we find that the relaxation obeys the simple aging in the 
asymptotic time regime. It is also very interesting to observe that the relaxation in the 
intermediate waiting times exhibit a sub-aging behavior (the characteristic relaxation time 
grows slower than the waiting time), which is often observed f30|in the experimental data of 
thermoremanent magnetization (TRM) of real spin glass systems. Note that the crossover 
from the sub-aging to the simple aging occurs over a rather broad time region. We find that 
close to the dynamic transition, as shown in Fig. 3(a), the simple scaling is not yet fulfilled 
for waiting times which are long enough to reach simple aging at T = 0.5. This suggests 
that the crossover time from sub-aging to simple aging becomes broader as the transition is 
approached. 

Though the sub-aging eventually makes a crossover to the simple aging behavior, due to 
the broadness of the crossover regime (particularly near the transition), we wanted to know 
whether there is a scaling function h(t), which can collapse the data in the intermediate 
times. One available empirical form is hit) = exp[j^(t/t ) 1 ^' 1 ], which has been successfully 
used to collapse the TRM data in spin glasses fl30f . The time t is the microscopic time 
scale (t = 1 for our case). The case of \i = corresponds to the absence of aging (no 
waiting time dependence), and the case of fi = 1 to the simple aging. Thus, the sub-aging 
behavior will give < fi < 1. As shown in Fig. 3(b), with this empirical form of hit) with 
fi ~ 0.81 for T = 0.61, we can collapse both the correlation and the integrated response with 
largest waiting times. We find that the exponent \x becomes larger in order to collapse the 
sub-aging data (with the same range of waiting times) at lower temperatures. We should 
emphasize that fi is an effective exponent, i.e., it tends to increase in order to incorporate 
the data with longer waiting times into scaling collapse due to the crossover to simple aging. 

In summary, we have developed an efficient method of integrating the dynamic equa- 
tions of a mean field spin glass model with p-spin interaction. This method allows us to 
see numerically for the first time the entire dynamic regime covering from microscopic to 
asymptotic regime. Tio cover such a broad dynamic range will also be of major importance 
in the sudy of aging in structural glasses f31| . We observe several new dynamic features 
in the intermediate time regime including the existence of a new time scale at which the 
breakdown of FDT sets in, the broad sub-aging regime and its crossover dynamics to the 
simple aging. These findings offer a new theoretical challenge for further understanding of 
the off-equilibrium dynamics of the p-spin model. 
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FIG. 1. (a) The correlation function C(t,t — r) and the integrated response function F(t, t — t) 
versus r with different t for T = 0.5. (b)The parametric plot —TF(t, t — t) versus C(t,t — r) with 
different t as in (a). The dashed straight line with slope —0.565 is shown for comparison. Inset: r p 
versus t in a logdog plot. The solid straight line with slope 0.68 is shown for comparison. 




FIG. 2. (a) The correlation function C(t w + r, t w ) versus r for different waiting times t w at 
T = 0.5. The upper dotted line is the numerical solution of (7). (b) C(t w + T,t w ) versus rescaled 
time t jt w for the data shown in (a). 
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FIG. 3. (&)C(t w +T, t w ) versus rescaled time r/t w for T = 0.61. (b) Scaling plot for C(t w +r, t w ) 
and F(t w + r, t w ) for T = 0.61 using h(t) = exp[t 1 -^ t /(l - p)] with n = 0.81. 
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